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The notion of an apparent horizon (AH) in a collapsing object can be carried over from the 
Lemaitre - Tolman (L-T) to the quasi-spherical Szekeres models in three ways: 1. Literally by 
the definition - the AH is the envelope of the region, in which every bundle of null geodesies has 
negative expansion scalar. 2. As the locus, at which null lines that are as nearly radial as possible 
are turned toward decreasing areal radius R. These lines are in general nongeodesic. The name 
"absolute apparent horizon" (AAH) is proposed for this locus. 3. As the envelope of a region, where 
null geodesies are turned toward decreasing R. The name "light collapse region" (LCR) is proposed 
for this region (which is 3-dimensional in every space of constant t); its boundary coincides with the 
AAH. The AH and AAH coincide in the L-T models. In the quasi-spherical Szekeres models, the 
AH is different from (but not disjoint with) the AAH. Properties of AAH and LCR are investigated, 
and the relations between the AAH and the AH are illustrated with diagrams using an explicit 
example of a Szekeres metric. It turns out that an observer who is already within the AH is, for 
some time, not yet within the AAH. Nevertheless, no light signal can be sent inside out through the 
AH. The analogue of AAH for massive particles is also considered. 
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I. MOTIVATION 

This paper deals with the relationship between ana- 
logues of an apparent horizon (AH) that exist in the 
quasi-spherical Szekeres models [l| - [2l| of the /?' ^ 
family. The AH was first defined by Szekeres @ as the en- 
velope of a region in which every bundle of null geodesies 
is converging. Hellaby and Krasinski gave the name 
of AH to a different entity, for which the name 'abso- 
lute apparent horizon' (AAH) is proposed here. The 
AAH is defined in terms of nongeodesic null lines that 
are, in a sense to be defined in Sec. [TV] as nearly ra- 
dial as possible. (Strictly radial curves do not exist in 
general Szekeres models because of their lack of sym- 
metry.) The AAH is the locus at which these nearly 
radial null curves are turned toward decreasing areal ra- 
dius R. In the Lemaitre ^ - Tolman [H (L-T) model, 
which is the spherically symmetric limit of the Szekeres 
models considered here, the AAH coincides with the AH 
24j , and the curves defining the AAH become radial null 
geodesies. 

One more analogue of AH results when we consider null 
geodesies and the region, where they are turned toward 
decreasing R. For this region, the name "light collapse 
region" (LCR) is proposed here. Unlike AAH and AH, 
which are 3-dimensional hypersurfaces in spacetime, the 



LCR is a 4-dimensional subset of spacetime because the 
family of geodesies defining it is not uniquely determined. 

The existence of the AAH is proven for every collapsing 
quasi-spherical Szekeres model in Sec. HVl In Sec, [V] the 
LCR is defined and it is shown that the 3-dimensional 
future boundary of LCR coincides with the AAH. In 
Sec. |VI] an explicit subcase of the quasi-spherical Szek- 
eres model is chosen for a detailed investigation. We 
illustrate the relation between the AH and the AAH 
by diagrams showing their positions in space. It turns 
out that, for some directions, an observer who is already 
within the AH is, for some time, not yet within the AAH. 
The analogue of AAH for massive particles is also con- 
sidered. In Sec. IVIII the matching of the quasi-spherical 
Szekeres solutions to the Schwarzschild solution is con- 
sidered. It is shown that it is the AH that matches to the 
Schwarzschild event horizon located at r — 2m and that 
the outgoing part of the AH is necessarily spacelike, so 
light rays cannot traverse it inside out. Both these facts 
indicate that the AH rather than the AAH is the true 
horizon. 

The aim of this paper is to gain more insight into the 
geometry of the Szekeres solutions. 



II. INTRODUCING THE /?' ^ 
QUASI-SPHERICAL SZEKERES SOLUTIONS 
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In this section, basic facts about the f3' ^ quasi- 
spherical Szekeres solutions are recalled for reference, and 
to define the notation. We will use the parametrization 
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introduced by Hellaby [25j . The metric of these solutions 
is 



F — f — 

t ~ 2 
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(2.1) 
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ds^ = dr 
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where E(z), P(z), Q{z), S(z) are arbitrary functions, and 
R(t, z) obeys the following equation; a consequence of 
Einstein's equations with dust source: 



R, t 2 = 2E(z) 



2M(z) 
R 



(2.2) 



M(z) being one more arbitrary function, and A being 
the cosmological constant. The coordinates of (|2.1I) are 
comoving, so the velocity field of the dust is = S^q 
and 11^=0. In the following we assume A = 0. 

This solution has in general no symmetry, and reduces 
to the L-T solution when P, Q and S are all constant. 1 
The sign of E(z) determines the type of evolution; with 
E(zq) < the matter shell at z = zq expands away from 
an initial singularity and then recollapses to a final singu- 
larity, with E(zq) > the shell is ever-expanding or ever- 
collapsing, depending on the initial conditions; E(zq) = 
is the intermediate case, ever-expanding with asymptot- 
ically zero expansion velocity, or its time-reverse. All 
three evolution types may exist in different regions of 
the same spacetime. We consider here the recollapsing 
(E < 0) solution of with A = 0, 
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where t_e(z) is one more arbitrary function and r/(t, z) is 
a parameter. The mass-density in energy units is 



2 (M, z -3M£, Z /£) 
R? (R, z -R£, z /£) ' 



(2.4) 



For e > 0, (M, z -3M£, Z /£) and (R, z -R£, z /£) must 
have the same sign. Note that the sign may be flipped 
by the transformation z — > —z, so we may assume that 

C = R lZ —R£, z /£ > at least somewhere. Let us then 
consider whether C can change sign as a function of z. 
The set where C = is either (1) a curvature singu- 
larity (a shell crossing - see a comment on it in the 



1 The "0' 0" refers to the fact that e* 3 d = R[£ depends on z, so 
0' = 0, z ^ (this notation follows Szekeres There exists 

another large family of Szekeres solutions, in which 0' = 0. They 
require separate treatment and will not be considered here. See 
an extended presentation in Ref. [lSl |. also for the associated 
quasi-plane and quasi- hyperbolic Szekeres models. 



next section) or (2), if it coincides with the set where 
M )Z —3M£, Z /£ = 0, an analogue of a neck, well-known 
from the L-T geometry [Isj . Case (1) is excluded by as- 
sumption - we choose the functions in the model so that 
shell crossings do not occur. In case (2), the neck (if it 
exists) is a global feature of spacetime, and no ordinary 
astronomical object in our neighbourhood is large enough 
to extend up to and through it. Thus, we assume that 
we are on one side of the neck, where 



R, z -R£ 1Z /£>0 
We also assume 



M, z -3M£, z /£ >0. (2.5) 



M, z >0 (2.6) 



because the region where M, z < occurs is an analogue 
of the region behind the equator of a closed space in the 
positive-curvature Robertson- Walker spacetimes. Again, 
no astronomical object is that large. 

The Robertson- Walker limit follows when z = r, 
R(t,z) = rS(t), E = Eor 2 , where Eq = constant and 
P = Q = 0, S = 1. This definition includes the defini- 
tion of the limiting radial coordinate (the Szekeres model 
is covariant with the transformations z = f(z'), where 
f(z') is an arbitrary function). 

The quasi-spherical model may be imagined as such a 
generalisation of the L-T model in which the spheres of 
constant mass were made non-concentric. The functions 
P{z), Q(z) and S(z) determine how the center of a sphere 
changes its position in a space t = constant when the 
radius of the sphere is increased or decreased [ll|. 

Within each single {t = constant, z = constant} sur- 
face, which is a sphere, the (x, y) coordinates of (|2.ip can 
be transformed to the spherical ($, tp) coordinates by 



(x-P,y- Q)/S = cot(tf/2)(cos</?,sin</>). 



(2.7) 



This transformation is called a stereographic projection. 
For its geometric interpretation see Refs. [l5| and [l~8l |. 
Using this transformation the factor £ , z /£ becomes 

£, z /£ = — [S, z cos-d + siwi)(P, z cosip + Q, z simp)] /S. 

(2.8) 



III. PROPERTIES OF THE QUASI-SPHERICAL 
SZEKERES SOLUTIONS 



Definitions of the Szekeres solutions based on invariant 
properties can be found using Ref. [l8j . 

Rotation and acceleration of the dust source are zero, 
the expansion is nonzero, the shear tensor is 



S = 



isdiag (0,2,-1,-1), 

R,tz —R,t R,z I R 
R, z —R£, z /£ 



where 



(3.1) 



The instant t = tsiz) in (|2.3[) is the Big Bang singu- 
larity corresponding to R — 0. When ts, z ^ (that is, in 
general) the instant of singularity is position-dependent. 



3 



Another singularity may occur where R, z —R£, z /£ = 
(if this equation has solutions for (x,y)). This is a 
shell crossing, but it is qualitatively different from that in 
the L-T model. As can be seen from (|2.1I) . the equation 
R, z —R£, z /£ — can define at most a subset of an {x, y} 
sphere. When a shell crossing exists, its intersection with 
a t = constant space will be a circle, or, in exceptional 
cases, a single point (in L-T it is a whole sphere). For 
more on shell crossings in all the Szekeres solutions see 
Refs. [lU and [2(|. They can be avoided if the functions 
and their derivatives obey a set of inequalities [H, [2l[ . 

Equation (|2.2p is formally identical with the Friedmann 
equation, but, with E and M depending on z, each sur- 
face z = constant evolves independently of the others. 

A quasi-spherical Szekeres region can be matched to 
the Schwarzschild solution across a z — constant hyper- 
surface 

The mass-density distribution given by (|2.4p can be 
decomposed into the spherically symmetric monopole 



2{M/ X f 



(R/ X ) 2 (R/xU 



where 



drf P 2 + Q 2 + S 2 + 1 

s 



and the dipole 

6MR, Z -2M, Z R X ,z-x£,z/£ 



R 2 (R, zX -Rx,z) R»-B£»/8' 



(3.2) 



(3.3) 



(3.4) 



The dipole is uniquely defined by the requirement that 
the surface where Ae = OJsure to exist, as follows from 
calculations - see Refs. [l6| and [Hj]) passes through the 
center of symmetry of the monopole. 2 



IV. APPARENT VS ABSOLUTE APPARENT 
HORIZONS IN THE QUASI-SPHERICAL 
SZEKERES MODELS 

The results of this section were partly reported in Ref. 
(l6j : the basic equations were introduced in Ref. [15j . 

An AH is the boundary of the region of trapped sur- 
faces. A trapped surface S is one on which the families of 
outgoing null geodesies on both sides of S converge (i.e. 
have a negative expansion scalar). Thus, if is any field 
of vectors tangent to null geodesies that intersect S, then 



s. 



Consequently, on an AH: 



(4.1) 
(4.2) 



2 Equation 1 13. 4t corrects a typo in eqs. (2.194) and (2.196) of Ref. 
[l6j and in eq. (19.165) of Ref. [lSj : one of the two appearances 
of <I> 2 in each equation should not be there. 



Proceeding from this definition, Szekeres Q found that 
in a quasi-spherical model the AH is given by the same 
equation as in the L-T model: 



R = 2M. 



(4.3) 



In an L-T model, (|4.2j) is equivalent to another definition: 
on an AH in collapsing matter, dR/dz calculated along 
an outgoing radial null geodesic changes from increasing 
to decreasing [18 1. 

Hellaby and Krasinski 15] considered the analogue of 
an AH in a quasi-spherical Szekeres model, using this 
second definition, but for non-geodesic null fields defined 
below. We propose to name this the "absolute apparent 
horizon" (AAH) - because even a maximally accelerated 
ray cannot get out of it. 

The reasoning was as follows. A general null direction 
k a = dx a /dt in the metric (|2.ip obeys 



(R, z -RE , z /£)' 



1 + 2E 



dzY 
dt) 



1 



£ 2 



dx 
~dt 



dy 
dt, 
'(4.4) 

Thus, on a null curve with dx/dt = = dy/dt (which, 
in general, will not be a geodesic [H, |dz/cl£| is 

maximal. Equation (14.4[) implies, along this path: 



df 

dz 



J 



VTT2E 



R,z 



R£,z 
~£~ 



j = ±l, (4.5) 



where j = +1 for outgoing rays, and j = —1 for ingoing 
rays. Intuition suggests that along a curve (|4.5I) the light 
signal should escape farther from the 'origin' R = than 
along any other path. An example in Sec. I VII will show 
that this is only partly true. Along some directions, the 
rays given by (|4.5p can indeed proceed toward larger val- 
ues of R even within the AH. But along the other direc- 
tions the reverse happens: the rays (|4.5p are re-directed 
to decreasing values of R even outside the AH. The rea- 
son for this behaviour is the fact that the two definitions 
of AH that are equivalent in the L-T limit are inequiv- 
alent in a Szekeres model: the locus where all bundles 
of light rays begin to converge is different from the locus 
where light rays are forced to collapse toward decreasing 
R, see Sec. El 

Let the solution of (14.51) be 



(4.6) 



t = t n (z). 



del' 



The value of R along this ray, R n (z) = R(t n (z), z), is 
a monotonic function of z in some neighbourhood of the 
emission point. The AAH is where R n (z) changes from 
increasing to decreasing or vice versa: 







dR n 
dz 



dR dt n dR 
dt dz dz 



J2M/W+2E( R£ 

75 . I R,z~ 



VT+2E 



£ 



R, 



(4.7) 



from (|2.2p and (|4.5p . where £ = +1 for an expanding 
model and I = — 1 for a collapsing model. We consider 
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an AAH that is created in the collapse phase (£ 
so it is defined by the outgoing rays (j = +1). 
(14.71) becomes 



R. 



( VT 



2E 



\ ^2M/R + 2E 



R 



£ J1 
£ 



0. 



= -1), 
Then, 



(4.8) 



In Ref. |15| it was found that in a constant-t space 
the AAH "is a kind of oval with half inside R — 2M and 
half outside" . 3 Here, we will investigate the relation of 
the AAH to the AH in a simple example of a recollapsing 
Szekeres model, and this will provide an illustration to 
the quoted statement. 

We use the following expression for R, z (to be calcu- 
lated from (EHJ); see eq. (18.107) in Ref. [lj]): 

3 This can be easily seen from 1)4. 81 . Suppose, for definiteness, 
that R, z > 0. Recall that £ > (evident from J27TJ). Then 



(S,z > 0) 



(R < 2M), (£, z = 0) => (R = 2M) and {£, z < 



Ell = (Mil _ Ell) + _ Eli 

R \ M E ) \2 E M 

sin 77(77 — sin?]) (— 2E) 3 / 2 



(1 — cos r]) 2 



M 



-t, 



sin 77 



(1 — COST?) 2 



(4.9) 



We note from (pOj) and (pOj) that with tt < 77 < 2tt, 
where R,t < 0, we have 



V / 2M/R + 2E = -V^2E 



sin 77 
1 — cos T] 



(4.10) 



We substitute (|4T9|) and (|4~T0| in (|4T8)l . then multiply the 
result by (1 — cost?) 2 to avoid the infinite values at 77 — » 
and 77 — ► 27r, and obtain: 

0) => (R > 2M). For the proof that £, z changes sign on an 
(x, y) sphere see Ref. [lf|; £, z = is a large circle on that sphere. 



dcf 



M, z 
M 



Eii 

E 



(1 - cost?) 3 / 2 + 



3^ 
2 £ 



— cos 7 



y/— 2Esinr] 



VI — cos ?y 



sin 77(77 — sin 77) (-2E) 
y/1 — cos 77 M 



3/2 



sm?7 
-^1 — cos 77 



V-2E sin 77(1 - cos 77)-^ = 0. 



(4.11) 



This determines 7?(-M, cc, ?/)aah. Then, from (|2.3j) : 
M 



i(M, x, y) 



AAH 



-2£) 3 / 5 



(?7-sin77) A AH + is- (4.12) 



We assume that shell crossings are absent. Among the 
conditions for no shell crossings, found in Ref. [l5j|, the 
following are useful here: 



2tt 



3 E, z M, z 



-2E) 



3/2 



2 E M J M 
(see eq. (126) in Ref. QJ]), and 

M, z /M - E, z /E>0 



t B . z < (4.13) 



(4.14) 



which follows from the fact that R, z /R > must hold 
for all (77,2), via (|4.9[) taken at 77 = n [l5j]. We observe 
that 



lim 



sin 77 



V^=lim^£=?==. (4.15) 



»7->-27r yjl — COS 77 ' »7^0 ^/l — COS 77 

Now we verify using (|4.11l) that 



being positive in consequence of (|4.14|) : and 
lim ^(77) 



(4.17) 



= 2V^2E 



. \^E- 

2 E 



M, z 
M 



(-2E) 



3/2 



<0, 



being negative in consequence of (|4.13l) . 

Thus *(7r) > and ^(2ir) < 0, so there exists an 
7/0 G (tt,2tt) at which ^(t]q) = 0, and it is unique (see 
Appendix [XJ. In passing, we have proved that each parti- 
cle in a re-collapsing quasi-spherical Szekeres model must 
cross the AAH before it hits the Big Crunch at 77 = 2tt. 



V. THE LIGHT COLLAPSE REGION (LCR) 
AND ITS FUTURE BOUNDARY 

Consider a bundle of geodesic light rays flashed simul- 
taneously from a common origin. Let v be the affine 
parameter along these rays, 9 be the expansion scalar 
of the bundle, k* 1 be the tangent vector field to the rays 
and 5S be the surface area of the propagating front of the 
bundle. Then the following holds f[lq. eq. (16.131)): 



=20 = — hx{5S). 
av 



(5.1) 



Consequently, by (|4.2[) . on an AH the SS stops increasing 
along the geodesies in the bundle and begins to decrease. 
In an L-T model, a light front (LF) flashed from the ori- 
gin R = remains spherically symmetric at all times, 
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and its surface area is proportional to R 2 calculated at 
the LF. Therefore, in a collapsing L-T model, the AH 
is at the same time the locus at which R\ LF reaches its 
maximum. This coincidence between 9 — and the max- 
imum of R\ L -p does not hold in a Szekeres model, as is 
demonstrated below. The LF is not spherically symmet- 
ric, different points on it have different R values at a 
constant t, so the area of the front is no longer propor- 
tional to R 2 . Thus, there may be locations where 9 < 0, 
but R is still increasing along the rays, and locations 
where 9 > while R is decreasing. This remark should 
help in understanding the relation between AH and AAH 
discussed in Sec. IVT1 

This non-coincidence allows us to define one more en- 
tity related to AH, which we propose to name the "light 
collapse region" (LCR). This is the region where R has 
extrema along null geodesies, and so, during collapse, the 
rays are turned toward the Big Crunch. 

Consider (|4.4[) along a null geodesic, and suppose we 
know the solution of the geodesic equations. Then (|4.4[) 
together with the geodesic equations defines the function 



t t n g( K Z') : 



(5.2) 



where "ng" stands for "along a null geodesic". For an 
outward-directed null geodesic we then have 



df 

dz 



R,z —R£,z /£ 



VT+2E U 



where 



1 - 



R? 
£ 2 



dx 
~d7 



(5.3) 



(5.4) 



>'<) 



Proceeding exactly as from (|4.6[) to (|4.8I) . and assum- 
ing we consider outward-directed null geodesies in the 
collapse phase of the model we arrive at the following 
analogue of 



dP 

dz 



ng 



P. 



y/2M/R + 2E 

~JT+WTF 



2E U 




\ V / 2M/R + 2E 
We define the LCR as the region where 



di? 



= 0. 



(5.5) 



(5.6) 



As stated before, the LCR is a 4-dimensional subset of 
spacetime and a 3-dimensional subset of a space of con- 
stant time. This is because the geodesies that define the 
LCR are not uniquely determined: the U in (|5.4[) depends 
on the direction of the geodesic considered and takes a 
range of values at a given (t, r). 

Consider the intersection of LCR with the AAH, i.e. a 
locus where (14.81) and (15.61) hold simultaneously. Within 



this set we have U = I => dx/dt = dy/dt = 0. Thus, 
LCR P| AAH is a set in which both the nongeodesic null 
curves referred to in (|4.5j) and the null geodesies referred 
to in (|5.3p - (|5.4[) begin to proceed toward decreasing 
R, and in addition the null geodesies happen to have 
dx/dt = dy/dt = there, i.e. to be tangent to the curves 
defining the AAH. The definition of this set is identical 
to the definition of AAH, eq. (|4.8j) . This shows that the 
AAH is the future boundary of the LCR. 



VI. EXPLICIT EXAMPLES OF AAH IN 
SZEKERES MODELS 

A. Null rays 

As an illustration, we take the recollapsing Szekeres 
model defined by the same equations that were used in 
Ref. [24[ to discuss the formation of galactic-size black 
holes in the L-T model: 4 



t B (M) = -bM 2 + t BQ , 
t c {M) = aM 3 +T Q + t B0 , 



(6.1) 
(6.2) 



where t B (M) is the bang time, tc(M) is the crunch time, 
a, b, t B Q and To are arbitrary constants; t B o is the time- 
coordinate of the central point of the Big Bang and To is 
the time between the Big Bang and Big Crunch measured 
along the central line M — 0. Then, from (|2.3p . since 
i] = 2?r at t = tc- 



2E(M) 



(2wM) 2 / 3 



(aM 3 + bM 2 + T ) 



2/3 ' 



(6.3) 



As shown in Ref. [15], eq. (185), the extreme values 
of £, z /£ are 



dot fj, 

£ 



± 



p» 



Q, 



s 



(6.4) 



However, in choosing (P, Q, S) precaution must be 
taken not to make D e too large. If it is too large, then ei- 
ther the numerator or the denominator of (|2.4[) becomes 
negative in a region of space, thus rendering the mass 
density negative there (and infinite where the denomi- 
nator is zero). Physically, this means that the dipole 
component of ()3.4p dominates over the monopole in that 
part of the space. 

We thus first choose such a value of D e that will make 
the difference between the graphs of AH and of AAH 
visible at the scale of a figure, and then we choose such 
P, Q and S that will imply the chosen value of D e . To 



4 The values of the parameaters a, b, To and tso used here will be 
different from those in Ref. l24ll . This model is meant to be an 
illustration to various geometrical possibilities; it is not supposed 
to describe any real object in the Universe. 
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maximize D e , at least one of the derivatives P, z , Q, z , S, z 
has to be large. Experiments showed that the following 
functions will yield the desired result: 

a = 0.1, b = 5000, T = 12.5, t B0 = 0, 
S = M - 29 , P = 0.5M - 29 , Q = 0, (6.5) 

and z' = M(z) was chosen as the new z-coordinate. The 
resulting AAH is shown in Fig. [TJ The figure shows t(M) 
on the AAH for two points in the (x, y) surface: the one 
where the D e given in (|6.4j) is maximal (positive) and 
where it is minimal (negative). These two curves are 
compared with the ordinary AH and with the crunch 
time function tp(M). (This figure is a modification of 
Fig. 1 in Ref. [2J|.) See Appendix [Bl for the proof that 
all four curves indeed have a common origin at M = 0. 
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FIG. 1: A comparison of the future absolute apparent horizon 
(AAH) with the ordinary future apparent horizon (AH) in 
the model defined by (|6.ip and (|6.2|) . Curve AAH+ is the 
AAH along the direction where the contribution from £, z /£ 
is minimal and curve AAH— is the AAH along the direction 
where this contribution is maximal. The dashed-dotted line 
represents the big crunch singularity. Horizontal solid lines 
show the instants for which the next figures are drawn, these 
are: t = 4.0, 10.0, 11.5, 11.9772, 12.3 and 13.0. Note that we 
use M as the radial coordinate. 

Figure [1] shows that the contribution from £, z /£ can 
either increase or decrease the region where the acceler- 
ated rays are forced toward the Big Crunch, depending 
on the direction. In the direction where this contribution 
is minimal (i.e. £, z /£ is negative - curve AAH—), the 
AAH appears later than the ordinary AH, and the term 
£, z /£ causes that the accelerating ray can still proceed 
toward increasing R in a region where a geodesic bundle 
already converges. In the direction where D e is maximal 
(i.e. £, z /£ is positive - curve AAH+), the AAH appears 
earlier than the AH, and the term £, z /£ causes that the 
accelerating ray is turned inward where a geodesic bundle 
is still diverging. 

Note what this means physically. A non-geodesic light 
ray is one that is guided by mirrors or optical fibers. 



When the AAH has a smaller radius than the AH, the 
observer who has already fallen into the AH still has a 
chance to send a message, using non-geodesic rays, to 
observers occupying loci with larger R. However, the 
nongeodesic ray has no chance to escape from inside the 
AH and will be turned toward the Big Crunch as well, 
only somewhat later than the geodesic one - see Sec. 
IVII1 Even this is possible only in some of the directions; 
in other directions the AAH is outside the AH and no 
ray within the AH, geodesic or not, can proceed toward 
larger R. (See later in this paper - Figs. [5]and|n]illustrate 
this point more clearly.) 

Figure [5] shows a 3d graph of M on the AAH as a 
function of x and y, at the time instant t = 13.0 (compare 
Fig. [IJ. All values of x and y are admissible, and at 
every pair (x, y) there will be an M obeying (|4. 1 1|) . The 
graph shows at which points in the (x, y) plane the AAH 
has the largest radius (as measured by M), and where 
it has the smallest radius. The maximum of M is at 
the intersection of curve AAH+ in Fig. Q] with the line 
t = 13.0; the minimum of M is at the intersection of 
curve AAH— with the same line. Comparison with Fig. Q] 
shows that the values of M are indeed all in the expected 
range. Figure [3] shows the intersection of the AAH with 
the ordinary AH, which, at t = 13.0 is at the mass 

M d = M AH = 3.82860 (6.6) 

The coordinates (x, y) are not very intuitive. To better 
visualize the AAH let us first use the stereographic pro- 
jection (12.71) to transform (x, y) to the (i9, ip) coordinates, 
and then map the AAH into an abstract Euclidean space 
with the coordinates (£, ip, £). The second transformation 
has the following form 

£ = Maah(&,<p) sintfcos^, 
ip = Maah($, f) sin^sintp, 
C = M AAH (0,<p)co80. (6.7) 

We now use these coordinates to present the evolution of 
the AAH. As seen from (|2.8p . when Q tZ = the extreme 
values of £ >z /£ with respect to ip are when ip = and ip = 
7T, which, as follows from (|6.7I) . implies ip = 0. Therefore, 
Fig. @] presents the intersection of the AAH with the 
plane (£, ip = 0, Q at 6 different time instants. 

Note how the lack of spherical symmetry influences 
the situation. The "origin", where R = (£ = 0, C = 0), 
is inside the smallest contour in Fig. 2) The AAH+ 
first appears off the origin (inside the closed curve in the 
lower left part of the figure). As seen, at this instant, 
most rays will miss it. Then it increases in diameter and 
encroaches on the origin. At the instant corresponding 
to the lowest point of AAH— in Fig. [T] (see the purple 
curve in Fig. [H which presents the AAH just moments 
before this instant, t = 11.9772), the cross-section is still 
connected, but consists of two tangent rings, one inside 
the other. The point of tangency lies at the minimum of 
AAH—. From that moment on, the cross-section splits 
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FIG. 2: Graph of M(x,y) on the AAH in the space t — 13.0 (compare Fig. [TJ. More explanation in the text. 




into two disjoint contours, the smaller of which becomes 
progressively smaller with increasing t, and shrinks to a 
point at the instant corresponding to the minimum of the 
Big Crunch (inside the smallest ring in Fig. [4j. 

Figure |4] does not show the cross-sections of the ordi- 
nary AH because they would obscure the image. The AH 
first appears shortly after t = 10.0 and, at the moment of 



first appearance, would show in Fig. |4] as a single circle 
with the center at (£, £) = (0, 0) and radius slightly larger 
than M = 2. At later instants, the cross-section of the 
AH splits into two circles, with the centers at the same 
point. The smaller circle has its radius decreasing as t 
increases, and shrinks to a point at t = Tq + tso = 12.5 
(the smaller contour of the AAH shrinks to a point at the 
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FIG. 4: The AAH at 6 different time instants (shown in Fig. [T]as horizontal lines). The curves shown are located on the plane 
(£, ip = 0, £)■ This figure can be imagined as a view of Fig. [T]by an observer sitting high on the i-axis and looking down; with 
one spatial dimension added. The value of Maah at a point (£o,Co) is the distance between (£o,Co) an d the origin R = 0, 
which is inside the smallest contour in the upper right area at (£, () — (0, 0). 



same instant). The larger circle of the AH keeps increas- 
ing, and intersects the larger contour of the AAH (which 
is not a circle) at two points at every instant. 

Three-dimensional surface-plots of the AAH and AH 
at t = 11.5 and t = 13 are presented in Figs. [3] and 
[6l At t = 11.5 the Big Crunch singularity has not yet 
appeared, the AAH is still a connected surface and the 
AH consists of two disjoint parts - one at Mah — 0.783 
and the other at Mah = 3.199. Both parts are presented 
in Fig. [S] (the dotted surfaces) , and, as seen, each one has 
one side inside the AAH and one side outside it (the solid 
surface). Each part of the AH has the shape of a sphere, 
while at t = 11.5 the AAH has the shape of a ping pong 
ball depressed on one side. 

At t — 13.0 the singularity already exists at R = 0. 
Each of AH and AAH consists only of a single surface, 
which surrounds the singularity. As before, part of the 
AAH is outside the AH, while the other part is located 
inside the AH. 

Between t = 11.5 and t = 13.0 there is a period when 
each of AH and AAH is split into two disjoint parts. 
We do not provide an illustration for this configuration 
because it would be unreadable. One can imagine it as 
the object from Fig. [5] that contains a small-scale copy 
of itself inside. The small object inside does not intersect 



with the large one. 

Figures 0] - [5] demonstrate that the AH and AAH do 
not in fact reveal the whole truth about the future fate of 
the light rays. There exists a region between the M = 
axis and the AAH— in Fig. [TJ in which future-directed 
rays are only formally not yet in the black hole because 
they are able to proceed outwards. However, they can 
do so only for a short while. They have no way to avoid 
intersecting the AAH and the AH in near future, and so 
they are doomed to hit the Big Crunch. 



B. Time-like curves 
Equation (|4.4j) for time-like trajectories is 



-Relief (**\ 2 = 

1 + 2E \ dt J 



' s 2 



For a trajectory with constant x and y we write 



(R, 



-RE, Z /Ef fdz » 2 
1 + 2E \dt' 



(6.8) 



(6.9) 
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FIG. 5: Left: The M-coordinate on the AAH (solid surface) and AH (dotted surface) represented as a function of # and ip 
in the space t — 11.5. The value of M{§,ip) is the distance of a point on the surface shown from the point (0,0,0), in the 
direction specified by (i?, The axes in this picture are in an abstract Euclidean space with coordinates (£,ip,C) used only 
to embed the AAH [see transformation ()6.7|) ]; the (£, VO do not coincide with the (x,y) of Fig. [2] The (M, ip) are spherical 
polar coordinates in this space. The AH consists of two disjoint spheres - one at Mah = 0.783 and the other at Mah = 3.199. 
The origin (£ = 0, ij) — 0, ( — 0) is inside the smaller AAH surface. Right: Intersections of the AAH (solid line), and the inner 
and outter AH (dotted lines) with the plane (£, ip — 0,() (analogous to Fig. \$$. 
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FIG. 6: The analogue of Fig. [5] at the later instant t = 13.0, when the singularity already exists at R — 0. 



where V = yjl- (|f) 2 < 1. Then gU becomes , R \/ 2M / R + 2E c jz /£ = Q 



V 



(6.10) 



R» (^/T+2E - \ /2M /R + 2E ^j and g-jT) becomes 
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*(»?) 



M, z E, z \ . , 3 / 2 

(1 — COS 77) ' 



cos 77 - 



\/-2£ , sin 



3^ 
2 £ 

»7 



M 



sin 77(77 — sin 77) (— 2E 1 ) 3 / 2 sinr; 

^ — — z — 

\/l — COS 77 M ' ^/l — COS 77 



V-^/l — cos 77 



V^2E . 

— - — sin 77(1 - cos 77) — = 0. 



(6.11) 
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FIG. 7: Comparison of the AAH of a null ray (dashed lines) 
with the AAH of a particle moving with velocity V = 0.9 (solid 
lines). The dotted line represents the AH, and the dashed- 
dotted line represents the big crunch (BC) singularity. As in 
Fig. [1] the inner curve (closer to BC) is AAH—, the outer 
curve is the AAH+. 



Figure [7] is the analogue of Fig. [T]for a particle moving 
with velocity V = 0.9. 



VII. WHICH IS THE TRUE HORIZON 
AAH? 



AH OR 



To get insight into this question we recall that a 
quasi-spherical Szekeres spacetime can be matched to the 
Schwarzschild spacetime across an z = b = constant hy- 
persurface; this was first proved by Bonnor Q. We 
recapitulate the basic facts about this matching by the 
method of Ref. [H. 

To verify the matching, the Schwarzschild solution 
must first be transformed to the Lemaitre [22J - Novikov 
(27j coordinates, see Ref. jl8[ (section 14.12) for a 
derivation. In these coordinates, it has the form 



ds 2 = dt 2 



R,r* 



1 + 2E(r) 



dr 2 - R 2 (t, r) (dd 2 + sin 2 dd<p 2 ) , 

(7.1) 



where R(t, r) is determined by the equation 

2m 
~R~' 



R, t 2 = 2E(r) 



(7.2) 



m being the Schwarzschild mass and E(r) being an arbi- 
trary function. In this form, the Schwarzschild metric is 
the limit M, r — of an L-T model, and the limit of con- 
stant M, P, Q, S of a quasi-spherical Szekeres solution. 

Further, the coordinates used on a sphere of constant 
(t, r) in (|7.ip must be transformed to those used in (|2.ip . 
Suppose the matching is to be done at r — z = b = 
constant. Then the transformation is 



t9 = 2 arctan ■ 



y/[x-P(b)]* + [y-Q(b)]* 



S(b) 



ip = arctan 



y-Q(b) 



x - P{b) 



(7.3) 



where P(b),Q{b) and S{b) are the values of the (P, Q, S) 
from (I2.1[) at z = b. The transformed metric (|7.1[) is 

ds 2 = dt 2 _ grr^ dr 2 _ ( da ,2 + dy 2^ ^ (? 4) 



1 + 2E(r) 



where 

c dof S{b) 
£■1 — 



x - P(b) 



S(b) 



V - Q{b) 



S(b) 



(7.5) 



Now it can be easily verified that the matching condi- 
tions between (|7.4p - (|7.5p and (|2.1I) are fulfilled at any 
r — z = b = constant, provided that the E(r) of (|7.5p 
and the E(z) of (|2.ip have the same value at r = z = b, 
and that both R(t, b) are the same function of t. The 
latter condition implies 



M(b) = m, 



(7.6) 



where M is the function from (|2.2p and m is the 
Schwarzschild mass from (17.21) . 

To answer the question asked in the title of this section 
we need to verify whether the AH is spacelike or other- 
wise. For the L-T models, this analysis was done in Ref. 
[1H and repeated in Ref. [l8j], with the result that the 
ingoing part of the AH (around R = in Fig. [1} can be 
any, while the outgoing part can be spacelike and point- 
wise null, but never timelike. We use the same method 
here, adapted to the Szekeres geometry, assuming that 
(12311 and (HU) hold. 

From (g3J we find R, t dt + R, z dz = 2M, Z dz along the 
AH, so 



dt 

dz 



2M, 2 -R, 



AH 



R,t 



(7.7) 



AH 
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FIG. 8: Top: The light cone at a point P on the AH and its 
intersection with the plane tangent to the (t, r) surface at P. 
The intersection determines a null direction ND - the thick 
line. The arrow in the plane marks a hypothetical direction 
tangent to the AH. Bottom: The plane from the upper panel. 
The figure shows also the past light cone of P (ND2) and three 
hypothetical directions of the vector tangent to the AH at P. 
In the position AH1, the AH would be outgoing timelike, in 
AH2 - spacelike (the position marked in the top panel), and 
with AH3 - ingoing timelike. If the direction of AH coincides 
with ND or ND2, then the AH is null at P. The quantity B 
defined in (|7.10|) identifies the various possibilities. 



But in the collapse phase R, t - 
along AH R = 2M, so 



s/2M/R +2E, and 



dz 



R,z ~2M, Z 



AH 



VTT2E 



(7.8) 



R=2M 



The equation of the AH is independent of (x,y), so a 
vector tangent to the AH has only the t- and r- compo- 
nents. We consider an intersection of a (i, r) surface with 
the light cone at a point of the AH, see Fig. HI Along a 
null geodesic that is tangent to this surface at the vertex 
of the light cone (and so has dx/dt — dy/dt — 0) we have 



dt 



R, z -2M£, z /£ 



ng/AH 



V1 + 2E 



(7.9) 



ng/AH 



As the lower panel of Fig. [8] shows, the following quantity 
indicates whether the AH is spacelike or otherwise 



B = 



(dt/dz) 



AH 



R, z -2M,. 



(di/d2) ng/ AH R, z -2M£, z /£ 
. 2 (M, z -M£, z /£) 



R, z ~2M£, z /£ 



ng/AH 

(7.10) 



ng/AH 



Namely 



THE AH IS 


WHEN 


outgoing timelike 


B > 1 


outgoing null 


B = 1 


spacelike 


-KB <1 


ingoing null 


B = -1 


ingoing timelike 


B < -1 



From (|23)1 and (|2l))) we see that M, z -M£, z /£ > and 
R, z —2M£, Z /£ > are always fulfilled, so necesarily 
B < 1, i.e. the AH can never be outgoing timelike or 
null; its outgoing part is necessarily spacelike. The other 
three possibilities listed in the table are allowed. 

As Fig. [T] shows, even if part of the ingoing branch of 
the AH is timelike, a null curve crossing the AH inside out 
will be trapped in the funnel formed by the AAH around 
R = 0. Whether it later crosses the AAH or not, it will 
be forced to hit the Big Crunch within a finite segment 
of its affine parameter. Where the AH is spacelike, a null 
line that crossed it once cannot cross it again without 
being re-directed toward the past. This shows that in its 
outgoing part the AH acts as a black hole surface, even 
in that region where the AAH is inside it. The conclusion 
is that the true horizon is the AH rather than the AAH. 

This conclusion is strengthened by the following con- 
sideration. If a portion of the Szekeres manifold, of fi- 
nite spatial diameter, is matched to the Schwarzschild 
solution, then, from (|7.6[) . the AH matches to the 
Schwarzschild event horizon. Thus, no signal can escape 
to infinity if it was within the Szekeres AH while crossing 
the outer surface at z = b. The intersection of the AAH 
with the outer surface of the Szekeres ball leaves no trace 
in the Schwarzschild geometry. In particular, this hap- 
pens in that part of the Szekeres region, where the AAH 
is earlier than the AH (e.g. in the left half of Fig. [T]). 

The above considerations are illustrated in Fig. |9l 
which shows a family of the nearly radial nongeodesic 
null lines (NRNL) for the model discussed in Sec. IVI Al 
It is similar to Fig. [TJ but for clarity we only show one 
part of it, where the AH is below the AAH. Calculations 5 



At each point of intersection of the NRNL with the AH and with 
the AAH we numerically calculated and compared their slopes. 
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FIG. 9: A family of the nearly radial null rays defined by 
(|4.5|l (thin solid lines). The solid thick line represents the 
A AH— , the dashed thick line represents the AH, and the 
dashed-dotted line represents the big crunch (BC) singularity. 
To the right of the dashed vertical line (M = 1.048) the AH 
is spacelike, to the left it is ingoing timelike. 



"kind of oval" is shown here in Figs. [5] and [6] 

An observer who would fall into the region between 
the surfaces of AAH and AH (top part of Fig. [5] be- 
tween the solid surface and dotted surface) for a short 
while would have a chance to send a message some dis- 
tance outwards (i.e on a path with increasing R). This 
is because in a general Szekeres model the hypersurface 
at which the light rays begin to converge (the AH) does 
not coincide with the hypersurface at which all rays are 
forced to proceed toward decreasing areal radius R (the 
AAH). However, the signal cannot proceed far enough 
to escape the AH. Moreover, if the Szekeres spacetime is 
matched to the Schwarzschild spacetime, the AH finds its 
prolongation in the Schwarzschild event horizon. Conse- 
quently, it is the AH that acts as a true horizon. 

Whether the concept of AAH can be usefully applied 
to astrophysical considerations about galactic black holes 
remains to be seen. For this purpose, the current position 
of the AH in space, inside the galaxy chosen for observa- 
tion, would have to be precisely determined, which seems 
to be a rather remote possibility. 



showed that, depending on the direction and the value of 
M, both the AAH and AH can be spacelike or ingoing 
timelike. Figure |H] shows the behaviour in the direction 
of maximal contribution from £_ z /£. In this direction the 
AAH— is everywhere spacelike, while the AH is ingoing 
timelike for M < 1.048 (the vertical line in Fig. |9|) and 
spacelike for M > 1.048. In other directions parts of the 
AAH can be timelike. The location of the border be- 
tween the ingoing timelike and spacelike parts of the AH 
is direction-dependent, as follows from (I7.10[) . (This is so 
because the slope of the null cone generator is direction- 
dependent, as seen from (|7.9p ). In neither case can the 
AAH or the AH be outgoing timelike. Thus, even if in 
some regions radial null rays can propagate toward in- 
creasing i?, they cannot escape from inside the AH and 
eventually they cross the AAH. 



VIII. SUMMARY 

In order to gain a deeper insight into the quasi- 
spherical Szekeres geometries we have investigated the 
spatial relation between the apparent horizon (AH) as 
denned by Szekeres and the absolute apparent horizon 
(AAH). The concept of the AAH was first introduced by 
Hellaby and Krasifiski [l5j], but under the name of AH. 
In Ref. [HI], this spatial relation was investigated at a 
general level, at which it was not possible to give graph- 
ical examples. Such graphical examples are given here 
in the simple subcase of the Szekeres model defined by 
(I6.1[) - (|6.2[) and (|6.5p . The examples illustrate what was 
said in Ref. [l5j], that the AAH "is a kind of oval with 
half inside R = 2M and half outside" . The shape of this 



Appendix A: The proof that the solution of (|4.11|) is 
unique 



We prove here that (j4.11j) is fulfilled by only one value 
of r] £ (it, 2tt) for each set of values of (z, x, y). 

We begin by recalling the following: 

1. Equation (|4.11[) that determines rj(z,x,y) on the 
AAH is derived from (|4.7I) . and so all quantities in it are 
calculated along the nearly radial nongeodesic null line 
(NRNL) that obeys (|4.4|) with x and y being constant, 



i.e. 



(R, z -R£, z /£) 2 fdz\ 2 
dt J 



1 + 2E 
2. Note, from (|4TT2"j) . that 



1. 



(Al) 



dt M 
dT] ~ (-2£) 3 / 2 



(1- cos 17) > for 77 e (0, 2tt), (A2) 



so t(rj) is a monotonic function in this range, i.e. 

(Vi < Vj) =>■ (KVi) < t(Vj )) for ever y (x,y). 

Now suppose that (|4.1ip has more than one solution for 
1] at a given (z, x, y), and call the solutions (771, ... , rjk), 
with rji < ■ ■ ■ < rjk- Then (|A2[) implies that there would 
be k instants t\ < ■ ■ ■ < t&, at which the given NRNL 
would intersect the AAH, all the ti,i = 1, . . . , k corre- 
sponding to the same (z, x, y) in (|4.11[) . Our supposition 
thus implies that the inverse function to t(z,x,y) has 
the property z(fi,x, y) = z(t2,x,y). Since this function 
is continuous (even differentiable, see (|A1[1 ). this means 
that for some t £ (£1,^2) we have dz/dt\ t _j = 0. But 
from (14.41) we have 



dz _ V1 + 2E 
dt ~ R, z -R£, z /£' 



(A3) 



13 



This can be zero only where E = —1/2. This set is a 
neck [l8| - an analogue of the Kruskal - Szekeres worm- 
hole in the Schwarzschild solution, i.e. a special location 
in spacetime that may or may not exist, depending on 
whether E attains the value —1/2 anywhere. 

Thus, for every NRNL obeying (|4.4|) that does not tra- 
verse a neck, dz/dt ^ everywhere along it. This means 
that z(t\,x,y) = z(t<i,x,y) cannot happen, i.e. that 
(14. lip has only one solution for r\ at each given [z, x, y). 
□ 



full range 77 G [0,27r]. We take M as the z-coordinate 
and observe from (I2.1[) that with (P, Q, S) given by (|6.5p . 
£,m /£ will be finite at M — 0. Then, substituting (|6.1j) 
and in (|4~TTj) we get 

lim 4- (77) = lim (1 ~ C ° S7?)2 = (B3) 

(because 'J (77) = all along the AAH). This is possi- 
ble only when limM-s-o cos 77 = 1, which, in the collapse 
phase, means 



Appendix B: The limit of AAH at the center M = 

For the numerical calculation we need to know the 
value of the function 7j(M) AAH in (|4~T2"j) at M = 0. This 
has to be calculated exactly because numerical programs 
are unreliable in calculating limits. From (|6.1[) and (|6.3[) 
we have 



lim t B (M) = t BQ , 

M->0 



lim t;(M)aah = t B o 



(B1) 



T 

— lim (77- sin 77). (B2) 
2vr iU-V ' i) \ i 



In order to calculate lmiM->o V we use (|4.11l) . Equa- 
tion (|4.15j) shows that sinTy/^l — cos 77 is finite in the 



lim 77 = 2tt. 



(B4) 



Using this in (IB2|) we obtain 



Hm i(M) AAH = t B o + T = Hm t c (M), (B5) 
i.e. the AAH at M — coincides with the Big Crunch. 
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